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' Abstract 

We propose a simple and systematic way of accounting for the back reaction on the background field 
due to the pair creation in the four-dimensional scalar QED. This method is straightforwardly generalizable 
to the gravity backgrounds. In the case of QED with the instantly switched on constant electric field 
background we obtain a remarkably simple formula for its decay rate. 
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1 Introduction 



Back reaction on classical background fields due to the quantum pair creation is an important problem for the 
understanding of the black hole and de Sitter space thermodynamics. The situation with the de Sitter space is 
the most controversial [Tj-|13j. At any rate, the proper way of taking into account the back reaction should be 
a first step towards quantization of gravity. Apparently the moment when the back reaction becomes relevant 
for the gravity background coincides with the moment when the gravity becomes quantum. 

All the existing ways of taking into account back reaction are rather complex and do not allow for a clear 
I/"") ' study of the situation. For that reason it is obviously tempting to analyze the phenomenon in a simple setting 
of the four-dimensional scalar QED with the electric field background, where the situation is supposed to be 
clear at least qualitatively. Even in the latter case the standard ways of accounting for the back reaction are 
rather involved (see e.g. [14] and [IS])- The goal of this note is to propose a simple way to calculate the decay 
rate of the background electric field, which can be straightforwardly generalized to the gravity backgrounds. 

Let us formulate the problem here. We assume that somehow at the moment of time t = a constant 
(everywhere in space) electric field was created. Saying another way, we neglect the boundary effects. We 
assume that at t = there are no electrically charged particles present on top of the field. At this moment one 
turns on interactions and the pair creation begins |16j (see as well [18] for the further development and |19j-|26] 
for the more recent progress and generalizations). On general physical grounds one expects that the value of 
the field will be decreasing at least due to the work performed to create pairs and to accelerate them. Due to 
the symmetry of the problem the decay should proceed homogeneously is space. 

The standard way to take into account the back reaction, considering the background field as classical, is to 
calculate the electric current which is created and to take into account the field due to this current [14] , |15) . 
If the field is directed along z axis one has to consider the Heisenberg evolution of the J3 component of the 
electric current up to some moment of time t and average it over the initial state of the problem. This way one 
would find the value of the created current at the moment t. Because the initial state is unstable one has to use 
the Keldysh-Schwinger diagram technic rather than the Feynman one |14j . Substituting the resulting value of 
the current into the RHS of the Maxwell's equation, one can find the field produced by it, which reduces the 
value of the original background field [14] . 

Such a procedure is rather involved and becomes even harder when generalized to the gravity backgrounds. 
In this note we propose another obvious and equivalent to the above mentioned one method to solve the problem 
in question. We propose to consider the Heisenberg-Euler effective action (including the Schwinger's imaginary 
contribution) for the background electric field in the scalar QED [TS], which is switched on at t = 0. Obviously 
to obtain the decay rate of the background electric field one just has to solve the equations of motion following 
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from the effective action. Along this line we find an unexpectedly simple answer (see (fH|) below). Our idea is 
to obtain the decay rate at the initial stage of the process and to approve that it is small for the subcritical 
values of the background electric field and as well to see whether the rate is big for the overcritical values of 
the electric field. Our eventual goal is of cause to find the decay rate for the overcritical field, because it is such 
a value of the field which is interesting from the perspective of the generalization of our considerations to the 
gravitational backgrounds. 

2 Effective action and equations of motion at one loop 

We would like to solve the following equation: 



d )l F fll ,(x,t) = (in j„(x,t) in 



(1) 



in the scalar QED. On the RHS of this equation stands the current due to the created pairs, which is calculated 
as the average of the Heisenberg evolution of the current operator j M over the initial state \in) [HI [15]. One 
can straightforwardly convert this expression into the interaction picture, where in our case the interaction part 
of the Hamiltonian is just the interaction of scalars with the classical electromagnetic field. Due to spacial 
homogeneity of the problem in question this average does not depend on x. 

The initial state \in) in question is the exact free theory vacuum state for the scalars in the constant electric 
field background. If such a state was an eigen-state of the Hamiltonian, i.e. if (out\S\in) , with S being the 
S-matrix operator, was a pure phase, then we could have made the following transformations: 
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where jy it) is the current operator in the interaction picture. Along these lines we would arrive at the 
standard functional integral in the in-out formalism. However, in our case the in-state is not an eigen-state of 
the Hamiltonian for all times and, hence, (out\S\in) is not a pure phase. Hence, we have to apply a different 
procedure due to Schwinger |16) . 

According to [17] we have to introduce the generating functional: 



>„. ■ < 3 » 

which involves a summation over a complete set of out -states. In our case, A± are gauge fields and, at the same 
time, sources for the currents (m|j M (i)|m) = 5W[A~, A + ]/SA+(t)\ A+=A - =A = -5W[A~, A+}/SA-(t)\ A+=A - =A . 
Note that the effective action obeys the conditions [T7]: W[A, A] = 0, W[A~ , A+] = -W[A+ , A ]*. 

The effective action in question for the background classical constant electric field in the four-dimensional 
scalar QED is equal to: 



s~ 



out, a ) ( out, a 



S 



log 



iW[A~, A+] 

+ 00 



[ V<f>*(x)V<t>(x) [ V<j>*_(x,t)V<f>-(x,t) expl- [ dt [ d 3 x (\D-<j)-\ 2 + m 2 \(t)-\ 2 )\ x 

xf Vct>* + {x,t)Vct> + {x,t) exp(~ / + dt fd 3 x (\D+cf> + \ 2 + m 2 \<j> + \ 2 )} 

JA,(t=0,S)=0 I 1 JO J V ' J 



(4) 



In this equation we have to impose (f>—(t = +oo,x) = = +oo,x) = 4>(x) at future infinity and assume the 
single common integration over this value of the field at the future infinity [T7] , <f>± are the scalar fields for the 
direct and reverse time directions, — 9 M — ieA^. 

With the use of the effective action under consideration, one can rewrite the equation ([T]) in the following 
way 
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d^F^ = i / V(f>*V(f>0*(ied„ + 2e 2 A„)0 expj-i f dt [ d 3 x 4>* (- D 2 + m 2 ) A = 

J <f>(t=0,g)=0 - 7 I ^ J J 



l<j>(t=0,x) = 

d d 
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In this equation we have unified </>± into one field <p an d according to (j4]) we have extended the time axis over 
the Keldish-Schwinger type contour C consisting of the two parts — C+ going from to +oo and C_ going 
backward. 

In the last line of ([5J, under the action of the differential operator, we have the in-in Keldish-Schwinger 
"++" type propagator on the half of R 31 . We would like to express this in-in propagator for the half of i? 3 ' 1 (t E 
[0, +oo)) through the in-out propagator for the full R 3 ' 1 (t G (— oo, +oo)). First, to express the in-in propagator, 
Gi n -in(y, z), through the in-out one, Gi n - ou t(y, z), for the full i? 3,1 (t € (— oo,+oo)), recall their definitions: 
G in -in(y,z) = (in\T<f>(y) (f>(z)\in) = 52 x f* n> x(y) fin,x( z ) and G in-out(y,z) = (out\T(f>(y) <f>(z)\in) / \out\in) = 

font \(v) fin.x(z) (out\b out .\ af n x \in) / \out\in) . In these expressions we have expanded the field <fr in terms of 
the in-harmonics fi n ,x and out-harmonics f ou t,\, which are properly normalized and correspond to the eigen- 
values A. Using Bogolyubov transformation from the out- to the in-harmonics one can express Gi n -i n (y, z) = 
Gin-out{y, z) + G* n _ out (y, z) — 2Re [Gi n - ou t{y, z)\. And finally, in order to express the in-out propagator for the 
full R 3 ' 1 space through the in-out propagator for the half-space t € [0, +00) with Dirichlet boundary conditions 
one just need to use mirror sources. 

To proceed we need to convert all the expressions into the Euclidian signature. We are going to make the 
Wick rotation back to the Minkowski space-time (xq = i t) at the end of the calculation. Using the considerations 
of the previous paragraph, it is straightforward to see that in Euclidian space the equation ([5]) can be rewritten 
as: 
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The point y is the position of the mirror source symmetric to y with respect to the surface xq = 0. 

The RHS of this equation can be calculated if the explicit expression for the background field is given. Hence, 
one has to look for the self consistent solution As(xo) of this equation. The solution should be self consistent 
in the sense that the RHS of ([6]), being calculated for the given ^(xo), should lead to the same ^(xo) as the 
solution of ([6]). Unfortunately, this problem has been solved only numerically [IB] . 

We would like to develop here the analytical approach in some approximation specified below. In particular, 
one can calculate exactly the RHS for the constant electric field background. The picture is self consistent if, as 
the result of ([6]) , the rate of the change of the background field is small. Otherwise one has to look for another 
approach. In any case our approximation works in the classical limit, when h — ¥ and the RHS of the equation 
in question is vanishingly small due to its one loop origin. 

So we assume that our background electric field is slowly changing in time and substitute into the LHS of 
([5]) F^v — — iE(xo)83\ fi 5„]o (imaginary unit is due to the Euclidean version of the theory). At the same time we 
neglect the time dependence of the field on the RHS, i.e. we put — ExqS^ into the propagator. Then the 
equation in question reduces to 
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In the case under consideration the RHS can be computed exactly. It is worth stressing at this point that the 
result of our calculation is gauge and Lorentz invariant because eq. ([6]) is written down in a gauge and Lorentz 
invariant way. 

We proceed with the computation of the propagator of charged scalars in the given external gauge potential: 
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z 2 ) 2 eE coth(eET) v , , s91 eE , w x , s 
— + ^ " [G/o - ^o) 2 + (y 3 - z 3 f] + i— (z + yo) (z 3 - y 3 ). (10) 
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is the classical action of the charged particle in the gauge field background in question. In the case of the 
propagator with y — » y one has to change only 

5cj [«,!/] = + 7 [(yo + zo) + [y 3 - z 3 ) \ +i— (z - y ) (z 3 - y 3 ). (11) 

Returning to the Minkowski space-time via the replacements £o — > it and E — > ii?, we get 

1 
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/n Tsin(e£ , T) i 

On general physical grounds one can expect only such a contribution to the RHS of (|12[) . which originates from 
the imaginary part of the Heisenberg-Euler Lagrangian describing the pair creation. 
Now, using that 
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we obtain 
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This formula gives us one loop exact answer for the decay rate of the background electric field in our approxi- 
mation. Ironically we find the exact decay rate of the field in the approximation when the field is held constant, 
which of cause is just a partial solution of the problem in any case. The generalization of our formulas to any 
dimension is straightforward. 

There are several points which are worth mentioning here. First, if we were considering the effective action in 
question in the full R 3,1 , as it was done in the original Schwinger's setting, rather than in the t > half of R 3 ' 1 , 
the RHS of (TTJ) would vanish, because of its time reversal (t — > —t) invariance. Second, it is straightforward to 
calculate that the VEV (Jo) in our setting does vanish as it should be, because the charges are created in pairs. 
Third, it is as well easy to observe that if one were considering the fate of the constant magnetic field instead of 
the electric one, he would have found that the field would remain constant, because the corresponding current 
does vanish. Which should be the case due to the absence of the imaginary contribution to the effective action 
in magnetic field background. 



3 Discussion and conclusions 

The decay rate (|14[) looks suspicious because of the dependence of its RHS on e 3 . One might think that as we 
change from particles to antiparticles the RHS changes the sign and, as the result, the background field grows 
rather than decays. We should stress that the RHS is always negative and depends on the modulus of the 
charge e. 
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To clarify the situation we would like to reobtain the answer (TUT) for the decay rate on the general physical 
grounds. Note that in our problem the initial density of charges in space is zero. Otherwise the space would 
have had non-zero conductance, a, and the decay rate would have been exponential: dE z /dt — j z — —aE z . 

Taking into account the set up of our problem, which is formulated at the end of the introduction section, 
we can see that the decay of the field proceeds homogeneously. Hence, the energy of the electric field per unit 
volume E 2 /8ir is spent on the work performed by the field on the creation of pairs. This work is proportional 
to e E z — e E t, taking into account that we set the light speed c = 1. Here z is the separation distance between 
the members of the pair reached during the observation time t. 

Let us clarify this relation. In the gauge when Aq is not zero, we have po — i m 2 + p 2 + p 2 + p 2 — eEz. 



Then, the distance between the turning points is Az — e ~^ x +Pv . Similarly in the gauge when A z is not 
zero, we have po — ^Jm 2 + p 2 + p 2 + (p z — eEt) 2 . Hence, the separation time between the turning points is 

At = — * — — g- 2 — - = Az. This observation establishes the relation between the separation distance and the 
observation time. 



Thus, 



^E 2 <x-2eEtw(E), (15) 



where w(E) oc e 2 E 2 e -1 ?^ is the approximate Schwinger's pair creation probability rate per unit time and unit 
volume. Thus, we restore the leading approximation of our exact one loop result (|14p. It is worth stressing at 
this point that somewhat similar formula (in the light cone frame) to (|15[) was obtained in |14) . 

The result (fl"4")l confirms the obvious expectation, that weak electric field (eE << m 2 ) changes slowly in time, 

because in this limit In ( 1 + e^ 32 ^ ) — > and therefore dE(t) jdt — > 0. Note as well that this result is applicable 



in case only when the decay rate is small, which is valid for the short period of time t << (m/eE) e~^~ . On 
the other hand, it gives us a hint about strong field limit, because naive (and illegal) application of this formula 
in the strong field limit (eE >> m?) tells us that there will be a fast decay E(t) oc 1/t 2 . I.e. in the case of the 
overcritical field we can not apply our approximation. 

How then one can find the back reaction rate in the case of the overcritical field? Should one still consider 
the background field as classical? Will the decay rate still be non-analytic in the value of the background field, 
if the field is grater than the critical value? The reason why we are interested in the overcritical electric field 
background is that it is the simpler model example predecessing the consideration of the gravity backgrounds, 
where the overcritical background field is the most interesting case. 
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